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Abstract 

Using recent measurements of the 6-quark fragmentation distribution obtained in 
e~^e~ — > bb events registered at the Z pole, the non-perturbative QCD component 
of the distribution has been extracted independently of any hadronic physics mod- 
elling. This distribution depends only on the way the perturbative QCD component 
has been defined. When the perturbative QCD component is taken from a parton 
shower Monte-Carlo, the non-perturbative QCD component is rather similar with 
those obtained from the Lund or Bowler models. When the perturbative QCD com- 
ponent is the result of an analytic NLL computation, the non-perturbative QCD 
component has to be extended in a non-physical region and thus cannot be described 
by any hadronic modelling. In the two examples used to characterize these two situ- 
ations, which are studied at present, it happens that the extracted non-perturbative 
QCD distribution has the same shape, being simply translated to higher-x values 
in the second approach, illustrating the ability of the analytic perturbative QCD 
approach to account for softer gluon radiation than with a parton shower generator. 



1 Introduction 



Improved determinations of the 6-quark fragmentation distribution have been obtained 
by ALPEH [1], DELPHI [2], OPAL [3] and SLD [4] collaborations which measured the 
fraction of the beam energy taken by a weakly decaying 6-hadron in e^e~ — >■ hb events 
registered at, or near, the Z pole. 

This distribution is generally viewed as resulting from three components: the primary 
interaction (e"'"e~ annihilation into a bb pair in the present study), a perturbative QCD 
description of gluon emission by the quarks and a non-perturbative QCD component 
which incorporates all mechanisms at work to bridge the gap between the prcvioTis phase 
and the production of weakly decaying 6-mesons. The perturbative QCD component can 
be obtained using analytic expressions or Monte Carlo generators. The non-perturbative 
QCD component is usually parametrized phenomenologically via a model. 

To compare with experimental results, one must fold both components to evaluate the 
expected x-dependence^ : 

t/ z z 

The final and the perturbative components are defined over the [0, 1] interval. As ex- 
plained, in the following, the non-perturbative distribution must be evaluated for x > 1, 
if the perturbative component is non-physical. The parameters of the model are then 
fitted by comparing the measured and predicted x-dependence of the 6-quark fragmenta- 
tion distribution. Such comparisons have already been made by the different experiments 
using, for the perturbative component, expectations from generators such as the JETSET 
or HERWIG parton shower Monte-Carlo. It has been shown, with present measurement 
accuracy, that most of existing models, for the non-perturbative part, are unable to give 
a reasonable fit to the data [1, 2, 3, 4]. Best results have been obtained with the Lund 
and Bowler models [11, 12]. 

In the following, a method is presented to extract the non-perturbative QCD com- 
ponent of the fragmentation function directly from data, independently of any hadronic 
model assumption. This distribution can then be compared with models to learn about 
the non-perturbative QCD transformation of 6-quarks into 6-hadrons. It can then be 
used in another environment than e'^e~ annihilation, as long as the same parameters and 
methods are taken for the evaluation of the perturbative QCD component. Consistency 
checks, on the matching between the measured and predicted ^-fragmentation distribu- 
tion, can be defined which provide information on the determination of the perturbative 
QCD component itself. 

In Section 2, the method used to extract the non-perturbative QCD component is 
presented. 

In Section 3, the extraction is performed for two determinations of the perturbative 
QCD component using: 

• the JETSET 7.3 generator [5], tuned on DELPHI data [6], running in the parton 
shower mode, 

1 A / ^ ^ X^ O 771 

^In the present analysis, x = ^ .-^ where Xe = —r- is the fraction of the beam energy taken by 

V rmn 

the weakly decaying hadron and Xmin = is its minimal value. 
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an analytic computation based on QCD at NLL order [7]. 



In Section 4 these results are discussed and insights obtained with the present analysis 
are explained. A parametrization for the non-perturbative QCD component is proposed. 



2 Extracting the aj-dependence of the non-perturbative 
QCD component 

The method is based on the use of the Mellin transformation which is appropriate when 
deahng with integral equations as given in (1). The Mellin transformation of the expres- 
sion for T>{x) is: 

ViN) = / dx x^-^ Vix) (2) 
Jo 

where A?" is a complex variable. For integer values of TV > 2, the values of V{N) correspond 
to the moments of the initial x distribution ^. For physical processes, x is restricted to 
be within the [0, 1] interval. The interest in using Mellin transformed expressions is that 
Equation (1) becomes a simple product: 

V{N) = Ppert.(A^) X Vnon-pertXN) (3) 

Having computed, in the A"-space, distributions of the measured and perturbative 

QCD components, the non-perturbative distribution, T) non-pert. {N) is obtained from 
Equation (3). Applying the inverse Mellin transformation on this distribution one gets 
T^ncm-pert.{x) without any need for a model input: 

I J AT '^■rneas.i^) -N 



T>non-pert.{x) = — f dN — — X (4) 

in which the integral runs over a contour in the complex A^-plane. The integration contour 
is taken as two symmetric straight half-lines, one in the upper half and the other in the 
lower half of the complex plane. The angle of the lines, relative to the real axis, is larger 
or smaller than 90 degrees for x values smaller or larger than unity, respectively. These 
hues are taken to originate from A^ = (1-01, 0). The contour is supposed to be closed by 
an arc situated at infinity in the negative and positive directions of the real axis, for the 
two cases respectively ^. 

In practice, the Mellin transformed distribution of present measurements, Vmeas.iN), 
has been obtained after having adjusted an analytic expression to the measured distri- 
bution in X, and by applying the Mellin transformation on this fitted function. The 
following expression, which depends on five parameters and gives a good description of 
the measurements (see Figure 1-left), has been used. 

D{x) =Po X \pixP^{l - xf^ + (1 -pi)xP^{l - xY^] (5) 



^By definition 25(1) (= 1) corresponds to the normalization of 'D{x). 

^It has been verified that tlic result is independent of a definite choice for the contour in terms of the 
slope of the lines and of the value of the arc radius. The result is also independent of the choice for the 
position of the origin of the lines, on the real axis, as long as the contour encloses the singularities of the 
expression to be integrated and stays away from the Landau pole present in 'Dpert. (■^) which is discussed 
in the following. 
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Figure 1: Left: Comparison between the measured (points with error bars) b- 
fragmentation distribution and the fitted parametrization using Equation (5). Right: 
Moments of the measured (full line) x distribution, of the perturbative QCD component 
[7] (dashed line) and of the generated distribution obtained in JETSET before hadroniza- 
tion (full line with circles). Data from [1] have been used. 



where po is a normalisation coefficient. Values of the parameters have been obtained by 
comparing, in each bin, the measured bin content with the integral of D{x) over the bin. 

Measurements of the 6- fragmentation distribution, in [1], have been published in a 
binned form, after unfolding of the experimental energy resolution. Values in the bins 
are correlated and, as the bin width is smaller than the resolution, the error matrix is 
singular. Only positive eigenvalues of this matrix have been used in the present analysis, 
when fitting parameters. 

The distribution of moments obtained with data from [1], and computed using the 
fitted distribution corresponding to Equation (5), is given in Figure 1-Right. The corre- 
sponding analytic expression is: 



P(iV)=po 



nP2 + N) , r(p4 + N) 

P\ — ; r + 1 — Pi — ; r 



(6) 



Quoted uncertainties, in Figure 1-Right, correspond to actual measurements and are 
highly correlated. They have been obtained by propagating uncertainties corresponding 
to the covariance matrix of the Pi,.,5 fitted parameters. When computing the moments, 
very similar results are obtained, for < 10, using directly the measured x-binned 
distribution. For higher values effects induced by the variation of the distribution 
within a bin, as expressed by Equation (5), have to be included. 

The Mellin transformed distribution of the JETSET perturbative QCD component 
has been obtained in a similar way, whereas the NLL QCD perturbative component is 
computed directly as a function of in [7] . At large values of A^, this last distribution is 
equal to zero for A^ = A"o — 41.7 and has a Landau pole situated at A^^, ~ 44. ^ 

"^Values for No and Nl depend on the exact values assumed for the other parameters entering into 



3 



3 a^-dependence measurement of the non-perturbative 
QCD component 

The X distribution of the non-perturbative QCD component extracted in this way depends 
on the measurements and also on the procedures adopted to compute the perturbative 
QCD component. In the following two approaches have been considered. The first one is 
generally adopted by experimentalists whereas the second is more frequent for theorists. 

3.1 The perturbative QCD component is provided by a gener- 
ator 

The JETSET 7.3 Monte Carlo generator, with values of the parameters tuned on DELPHI 
data registered at the Z pole has been used. Events have been produced using the parton 
shower option of the generator and the 6-quark energy is extracted, after radiation of 
gluons, just before calling the routines to create a 6-hadron that takes a fraction z ^ 
of the available string energy. This distribution is displayed in Figure 2. It has to be 
complemented by a 5-function at a; = 1 which contains ~ 4% of all events. In this peak, 
6-hadrons carry all the energy of the 6-quark as no gluon has been radiated. 

Applying the method explained in Section 2, the corresponding non-perturbative QCD 
component has been extracted, and is displayed also in Figure 2. Above x = 1, it is 
compatible with zero, as expected. The quoted error bar, for a given value of x, has 
been obtained by evaluating the values of non-pert. {^) for different shapes of the 6-quark 
fragmentation distribution which are obtained by varying parameters pi,.,5 according to 
their measured error matrix. 

3.2 The perturbative QCD component is obtained by an ana- 
lytic computation based on QCD 

The perturbative QCD fragmentation function is evaluated according to the approach 
presented in [7]. This next to leading log (NLL) accuracy calculation for the inclusive 
6-quark production cross section in e^e" annihilation, generalises previous calculations 
by resumming the contribution from soft gluon radiation (which plays an important role 
at large x) to all perturbative orders and to NLL accuracy. These computations are 
done directly in the A^-space. Soft gluon radiation contributes to the logarithm of the 
fragmentation function large logarithmic terms of the type (logA^)^, with p < n + 1. 
These terms appear at all perturbative orders n in a^. In the calculation at NLL accuracy 
[7], the two largest terms, corresponding to p = n + 1 and n, have been resummed at 
all perturbative orders. The calculation is expected to be reliable when A^ is not too 
large (typically less than 20). To obtain distributions for the variable x from results 
in moment space, one should apply the inverse Mellin transformation, that consists in 
integrating over a contour in A^ (Section 2). When x gets closer to 1, large values of A" 

the computation; see Section 3.2 where values of these parameters have been hsted. 

— ^b^pt" is the boost-invariant fraction of the 6-jet energy taken by the weakly decaying B meson. 

This variable is defined for a string stretched between the 6-quark and a gluon, an anti-quark or a diquark. 
In the present analysis no distinction has been introduced between z and x as no string model has been 
considered. 
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Figure 2: x- dependence of the perturbative (dotted line) and non-perturbative (full line) 
QCD components of the measured [1] b-fragmentation distribution. These curves are 
obtained by interpolating corresponding values determined at a large number of points 
in the x-variable. Quoted error bars correspond to measurement uncertainties and are 
correlated for different x-values. The perturbative QCD component is extracted from the 
JETSET 7.3 Monte Carlo generator. The dotted curve has to be complemented by a 
6-function containing 4% of the events, located at x = 1. 

contribute and thus the perturbative fragmentation distribution is not rehable in these 
regions. This behaviour affects also values of the distribution at lower x as moments of 
this distribution are fixed. In addition to the break-down of the theory for large values 
of A^, uncertainties attached to the determination of the theoretical perturbative QCD 
component are related to the definition of the scales entering into the computation. This 
component also depends on two parameters: the 6-quark pole mass (mi,) and Ag^^,, that 

have been taken as mf'^ = (4.75 ±0.25) GeV/c^ and Ag^^ = (0.226 ± 0.025) GeV. Scale 
and parameter depending variations of the moments of the perturbative QCD component 
are given in Figure 3. These variations are fully correlated versus A^. 

The extracted non-perturbative component is given in Figure 4. Its shape depends 
on the same quantities as those used to evaluate the perturbative distribution, and thus 
similar variations appear, as drawn also in the Figure. 

It has to be noted that the data description in terms of a product of two QCD compo- 
nents, perturbative and non-perturbative, is not directly affected by uncertainties attached 
to the determination of the perturbative component. This is because the non-perturbative 
component, as determined in the present approach, compensates for a given choice of 
method or of parameter values. 

To obtain the complete expected x-distribution of 6-hadrons from the theoretical cal- 
culation, one has to be able to evaluate the integral given in Equation (1). When x 
becomes close to 1 there are numerical problems, and consequently, the high-x (x > 0.96) 
behaviour of the perturbative QCD component has been studied. As mentioned before, 
this region corresponds to high-N values where the perturbative approach fails. As a 
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Figure 3: Variations of the calculated perturbative QCD component [7], depending on 
the renormalization scales (fi, fj,o), the factorisation scales (fip, fJ'Op), the b-quark mass 
(nib) o-nd AqIjj^. The full lines are corresponding to the central values (fJ^ = fJ^F = Q = 
91.2 GeV, /io = /ioF = mb, rub = 4.75 GeV/c^ and AgJ.^ = 0.226 GeV;. 

result, the high-x behaviour of the distribution is non-physical; it oscillates. To have 
a numerical control of the distribution in this region, it has been decided to take into 
account x values which are below a given maximum value, Xmax, above which the distri- 
bution is assumed to be equal to zero. Moments of this truncated distribution show a 
small discrepancy when compared with moments of the full distribution. This difference 
has an almost linear dependence with A^. To correct for this effect, x^ax is chosen such 
that the difference between moments is a constant value (the slope in N being close to 
zero at this point). This difference can then be corrected by adding simply a 5-function 
at a; = 1, so that the total distribution is normalized to 1. A typical value for x^ax is 
0.997 and the 6 component corresponds to 5% of the distribution. We stress that the 
truncation at Xmax and the added ^-function do not contribute in the determination of 
the non-perturbative component using Equation (4). But this procedure is necessary for 
checking that the extracted non-perturbative component in the x-space, when convoluted 
with the perturbative distribution, effectively reproduces the measurements and also for 
testing hadronic models given in the x-space. 

Conversely to the perturbative QCD component which was, by definition in [7], defined 
within the [0, 1] interval, the non-perturbative component has to be extended in the region 
X > 1. This "non-physical" behaviour comes from the zero of Pperi.(A^) for N = Nq which 
gives a pole in the expression to be integrated in Equation (4). Using properties of 
integrals in the complex plane, it can be shown that, for x > 1, the non-perturbative 
QCD distribution can be well approximated by x~^°. 

Errors bars, given in Figure 4, have been obtained using the same procedure as ex- 
plained in Section 3.2. 
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Figure 4: x-dependence of the perturbative (dotted line) and non-perturbative (full line) 
QCD components of the measured b-fragmentation distribution. These curves are obtained 
by interpolating corresponding values determined at a large number of points in the x- 
variable. The perturbative QCD component is given by the analytic computation of [7]. 
The thin lines on both sides of the non-perturbative distribution are corresponding to 
A^o = /^OF = {"^6/2,2mfc} (dotted lines) and Ag^^ = (0.226 ± 0.025) GeV (dashed lines). 
Variations induced by the other parameters, = fip = {Q/2,2Q} and mi, = (4.75 ± 
0.25) GeV/c^ are smaller. In addition, quoted error bars correspond to measurement 
uncertainties and are correlated for different x -values. The perturbative QCD dotted curve 
has to be complemented by a 6 -function containing 5% of the events, located at x = 1. 

4 Results interpretation 

The x-dependence of the non-perturbative QCD component, obtained in this way, does 
not depend on any non-perturbative QCD model assumption but its shape is tightly 
related to the procedures used to evaluate the perturbative component, and thus, the two 
distributions have to be used jointly. 

4.1 Comparison with models 

Non-perturbative components of the 6-quark fragmentation distribution, taken from mod- 
els, have been folded with a perturbative QCD component obtained from a Monte Carlo 
generator and compared with measurements by the different collaborations [1, 2, 3, 4]. 
Apart for the Lund and Bowler models, none of the other models provided a good fit to 
the data. In Figure 5 the directly extracted non-perturbative components are compared 
with distributions taken from models [8, 9, 10, 11, 12] whose parameters have been fitted 
on data from [1]. Results have been obtained by comparing, in each bin, the measured 
bin content with the integral, over the bin, of the folded expression for Vpredictedix) . They 
are summarized in Table 1. It must be noted that parameters given in this Table, when 
the perturbative QCD component is taken from JETSET, may differ from those quoted 
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in original publications as an analytic computation is done in the present study, using 
Equation (1), whereas a string model is used in the former; other sources of difference 
can originate from the exact values of the parameters used to run JETSET and from the 
definition of x which varies between and 1 in our case. Numerically, if one compares 
present results with those obtained by OPAL [3], in their own analysis, for the values of 
the parameters of the Lund and Bowler models ^, differences are well within uncertainties. 



Model 


JETSET 


NLL Pert. QCD 




par am. 


X'' /NDF 


par am. 


/NDF 


Kartvelishvili [8] 
x^b(l - x) 




27/4 


£(, = 14.0 ± 0.7 


22/4 


Peterson [9] 

X \ X i.—X J 


= {5-6+^;^) X 10-3 


62/4 


Eft = (3.1 ± 0.5) X 10-3 


29/4 


C.S [10] 


61, = p X 10-3 


144/4 


«6 = (4.4+;;^)xio-* 


18/4 


Lund |11J 


a = 1.55 ± 0.15 


2.6/3 


a = 0.04 ± 0.09 
^"^bl. = 8-8-0.8 


0.5/3 


Bowler [12] 

5--5— (1 - x)<»exp 1 !ii 1 


a = 0.77t-;^ 


2.5/3 


a = 0.00 ± 0.03 

hm^^ = 53. =t 6. 


4.4/3 



Table 1: Values of the parameters and of the 'x^ /NDF obtained when fitting results from 
Equation (1), obtained for different models of the non-perturbative QCD component, to the 
measured b-fragmentation distribution. The two situations corresponding, respectively, to 
the perturbative QCD component taken from JETSET or from [7] have been distinguished. 
The Lund and Bowler models have been simplified by assuming that the transverse mass 
of the b- quark, mb±, is a constant. 

When the perturbative QCD component is taken from the analytic NLL calculation 
the Lund and Bowler models give also the best description of the measurements. Values 
for the parameter a in these models arc even compatible with zero corresponding to a 
behaviour in l/z°' exp (—A/z), to accomodate the non-zero value of the non-perturbative 
QCD component at x = 1. 

But, as these models have no contribution above x — 1, their folding with the pertur- 
bative QCD component cannot compensate for the non-physical behaviour of the latter. 
The folded distribution is oscillating at large values (see Figure 6-Left). In particular, 
the predicted value in the last measured x-hin, which is found to be in reasonable agree- 
ment with the measurements after the fitting procedure when using the Lund and Bowler 
models, results from a large cancellation between a negative and a positive contribution 
within that bin. This is more clearly seen when considering the moments of the overall 
distribution, which are given by Equation (3). For moments of order A^, the weight x^~^ 
introduces a variation within the x-hin size, which was not accounted for in the previous 
fit in x, and effects are amplified mainly at large N values which correspond to the high 
X region. This is illustrated in Figure 6-Right. 

From this study, it results that all models have to be discarded when folded with 
the NLL perturbative QCD fragmentation distribution. The goodness of the fit in x, as 

®This numerical comparison has been restricted to these two models as other models do not give a 
good description of the measurements. In addition, only OPAL has fitted the parameters of the two 
models assuming that the transverse 6-mass was a constant, as done also in our analysis. 
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measured by the corresponding value, does not reflect all the information because it 
was not required that the folded distribution remains physical (positive) over the [0, 1] 
interval. This folding procedure has thus to be considered only as an exercise and the 
non-perturbative QCD distribution has to be extracted from data. 




1 .2 



Figure 5: Comparison between the directly extracted non-perturbative component (thick full 
line) and the model fits on data taken from [1]. Left: the perturbative QCD component is 
taken from JETSET. Right: the theoretical perturbative QCD component[7] is used. 




Figure 6: Comparison between the measured and fitted x -distributions using different mod- 
els. Left: the measured binned- distribution in x is compared with fitted results. Right: 
moments of the corresponding distributions are compared. 
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4.2 Proposal for a new parametrization 

As explained in Section 2, the non-perturbative QCD component of the 6-fragmentation 
distribution, has been extracted independently of any hadronic model assumption but it 
depends on the modelling of the perturbative QCD component. 

When a Monte-Carlo generator is used to obtain the perturbative component, it can 
be verified if the non-perturbative component has a physical behaviour for all x-values. 
In Figure 2, below x = 0.6, (4 ± 3)% of the integrated distribution is negative. A larger 
deviation would have indicated some inconsistancy between experimental measurements 
and gluon radiation, as implemented in the generator. 

Such a test cannot be made, a priori, when the perturbative QCD component is taken 
from an analytic computation as this distribution is already unphysical in some regions. 
The non-perturbative extracted distribution, as given in Figure 4, is precisely expected 
to compensate for these effects. It can be noted that, also in this case, the distribution is 
compatible with zero below x = 0.6. This shows that the perturbative QCD evaluation 
of hard gluon radiation is in agreement with the measurements. The small spike, close 
to X = 0, related to the multiplicity problem [13] in the perturbative evaluation, has no 
numerical effect in practice. 

To provide an analytic expression, which agrees better with the extracted point-to- 
point non-perturbative QCD distribution, the following function has been used: 

T^non-pert.{x) = -/V(p) exp ^ for X < Xq] 

= N{j>) exp— ^ for xq < x <1\ 

= A^(p)x-^° for X > 1; p = (xq, (7_, (7+) (7) 

It corresponds, for < a; < 1, to Gaussian distributions with different standard deviations 
when X is situated on either sides of xq. As explained in Section 2, the behaviour of the 
non-perturbative distribution for x > 1 is related to the presence of a zero in T>pert.{^) 
located at = A^^q. When the perturbative distribution is taken from a Monte-Carlo, 
there is not such a zero and it has been considered that 'Dnon-pert.{^) — when a; > 1. 
When the perturbative distribution is taken from theory, the value of Aq has been fixed 
to 41.7. The value for Aq depends a lot on the central values for the parameters and the 
scales, adopted in the perturbative QCD evaluation ^: 

A^o - 41.7;i^(Agi^) ± 2.5im,)+_lt'A^io = /^of) (8) 

A^(p) is a normalisation factor such that the integral of the expression, given in Equa- 
tion (7), between and oo is equal to unity. 

The parameter values, obtained when the perturbative QCD component is taken either 
from JETSET or from theory, are given in Table 2. It can be noted that the shape of 
the non-pcrturbativc QCD distribution is similar in the two cases, the maximum being 
displaced to an higher value in the latter. 

Comparison, in x- and in moment-space, between the measured and fitted distributions 
are given in Figure 7. For moments higher than 40, the fitted moments strongly deviate 
from the measurements as one is approaching the Landau pole and the formalism of [7] 
is no longer valid. 

''A^o is independent on the value assumed for the other scale fjb = fip- 
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^+ 


XVNDF 


JETSET 


0.911 ±0.006 


0.088 ±0.006 


0.045 ±0.005 


3.0/2 


NLL pert. QCD 


0.955 ±0.011 


0.087 ±0.006 


o.o6l[!:^^ 


0.7/2 



Table 2: Values for the fitted parameters of the non-perturbative QCD component corre- 
sponding to JETSET and NLL perturbative QCD using data from [1]. 
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Figure 7: Comparison between the fitted and measured b-quark fragmentation distribu- 
tions. Left: Differences between the fitted and the measured distributions in each x-bin 
are shown and the total error bars, relative to the measurements, are displayed. The fitted 
results correspond to the averaged values obtained from Equation (1) over the bin. Right: 
Differences between the fitted and measured moments of the h-fragmentation distribution. 
Fitted moments result from the product given in Equation (3) in which V non-pert. [N) are 
the moments of the fitted non-perturbative QCD distribution corresponding to Equation 
(7) and VpertX^) corresponds to the analytic computation [7]. 

5 Conclusions 

The measured 6-quark fragmentation distribution has been analysed in terms of its per- 
turbative and non-perturbative QCD components. 

The x-dependence of the fragmentation distribution has been extracted in a way which 
is independent of any model for non-perturbative hadronic physics. It depends closely on 
the way the perturbative QCD component has been evaluated. The obtained distribution 
differs markedly from those expected from various models. 

Below X = 0.6, this distribution is compatible with zero indicating that most of gluon 
radiation is well accounted by the perturbative QCD component evaluated using the 
LUND parton shower Monte-Carlo or computed analytically. 

As the non-perturbative QCD distribution is evaluated for any given value of the x- 
variable it can be verified if it remains physical over the interval [0, 1] when used with 
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a Monte-carlo generator which provides the perturbative component. The evidence for 
unphysical regions would indicate that the simulation or the measurements are incorrect. 
There is not such an evidence in the present analysis. 

Above X = 0.6, the obtained distribution is similar in shape with those expected from 
the Lund symmetric [11] or Bowler [12] models, when the perturbative QCD component 
is taken from JETSET. This is no longer true, also for these two models, when the 
perturbative QCD component is taken from the analytic result of [7] . It has been found 
that, because of the analytic behaviour of the perturbative QCD component, the non- 
perturbativc QCD distribution must be extended above x = 1. The x-behaviour of the 
non-perturbative component, for ,x > 1, is determined by the possible existence of a zero 
in T>pert.{^)j for A*" > 0. When the perturbative component has non-physical aspects, 
it is thus not justified to fold it with any given physical model. An approach has been 
proposed to solve this problem and a parametrization of the obtained distribution has 
been provided. 

The non-perturbative component, extracted in this way, is expected to be valid in a 
different environment than e~^e~ annihilation, as long as the perturbative QCD part is 
evaluated within the same framework (analytic QCD computation or a given Monte Carlo 
generator), and using the same values for the parameters entering into this evaluation as 
m^°'^, ^QCD or generator tuned quantities. It is planned to do a similar analysis using 
the improved analytic computation of the perturbative QCD 6-fragmentation distribution 
given in [14]. 
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